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Q \ Abstract. We study the behavior of the pion dispersion relation in a pion medium at finite density 

D ■ and temperature, introducing a chemical potential to describe the finite pion number density. Such 

' description is particularly important during the hadronic phase of a relativistic heavy-ion collision, 

between chemical and thermal freeze-out, where the pion number changing processes, driven by 
the strong interaction, can be considered to be frozen. We make use of an effective Lagrangian 
that explicitly respects chiral symmetry through the enforcement of the chiral Ward identities. The 
pion dispersion relation is computed through the computation of the pion self-energy in a non- 
^ ' perturbative fashion by giving an approximate solution to the Schwinger-Dyson equation for this 

\ self-energy. The dispersion relation is described in terms of a density and temperature dependent 

■ mass and an index of refraction which is also temperature, density as well as momentum dependent. 
CO ' The index of refraction is larger than unity for all values of the momentum for finite pL and T . Given 

, the strong coupling between p vectors and pions, we argue that the modification of the pion mass 

■ due to finite pion density effects has to be taken into account self-consistently for the description of 
I the in-medium modifications of p's. 

INTRODUCTION 

Oh! 

^ . In recent years, the possibility to produce a locally thermalized, deconfined state of mat- 

ter, whose degrees of freedom are the quarks and gluons of quantum chromodynamics 
(QCD), the so called quark-gluon plasma (QGP), in high-energy heavy-ion collisions, 
has attracted a great deal of attention, both on the experimental and the theoretical as- 
pects of the subject [|l]]. 

If the QGP is produced in these kind of reactions, the prevailing view portrays the 
evolution of such a system traversing a series of stages, the last of which consists of 
a large amount of hadrons strongly interacting in a finite volume, until a final freeze- 
out. For not too high temperatures, hadronic matter consists mainly of pions, therefore, 
the study of the propagation properties of pions within the above described environment 
represents an important ingredient for the understanding of the properties of the hadronic 
system at and just before freeze-out [@, Q]. 

The hadronic degrees of freedom are customarily accounted for by means of effec- 
tive chiral theories that incorporate the Goldstone boson nature of pions. One of such 
theories is chiral perturbation theory (ChPT) which has been employed to show the well 
known result that at leading perturbative order and at low momentum, the modification 
of the pion dispersion curve in a pion medium at finite temperature is just a constant, 
temperature dependent, increase of the pion mass [0]. ChPT has also been used in a 



two-loop computation of the pion self-energy ^ and decay constant [Q] . A striking re- 
sult obtained from such computations is that at second order, the shift in the temperature 
dependence of the pion mass is opposite in sign and about three times larger in mag- 
nitude than the first order shift, already at temperatures close to 150 MeV. This result 
might signal either the breakdown of the perturbative expansion at these temperatures or 
the need to perform such calculations using schemes other than perturbation theory. 

However, a missing ingredient in the calculations of the pion dispersion curve is 
the treatment of the medium's finite density. The conceptual difficulty is related to the 
fact that, though it is possible to assume that the system is in (at least local) thermal 
equilibrium, strictly speaking the only conserved charge that can be associated to the 
pion system is the electric charge and thus, for an electrically neutral pion system 
the corresponding chemical potential vanishes. The behavior of the pion mass in the 
presence of an isospin chemical potential has been recently studied in Ref. [Q]. But in 
order to describe a situation in which the number of pions in thermal equilibrium is 
finite, we need to consider a chemical potential associated with the pion number, instead 
of its charge. 

Recall that the pion number is not a conserved quantity due to either strong, weak 
or electromagnetic processes. Nevertheless, the characteristic time for electromagnetic 
and weak pion number-changing processes, is very large compared to the lifetime of 
the system created in relativistic heavy-ion collisions and therefore, these processes 
are of no relevance for the propagation properties of pions within the hadronic phase 
of the collision. As for the case of strong processes, it is by now accepted that they 
drive pion number toward chemical freeze-out at a temperature considerably higher 
than the thermal freeze-out temperature and therefore, that from chemical to thermal 
freeze-out, the pion system evolves with the pion abundance held fixed ^. Under 
these circumstances, it is possible to ascribe to the pion density a chemical potential and 
consider the pion number as conserved [ |TOt [TT] ]. In this context, the role of a finite pion 
chemical potential into a hadronic equation of state has been recently investigated in 
Ref. [0. 

Furthermore, another important ingredient in the analysis is the well know fact that in 
finite temperature field theories with either massless degrees of freedom or that exhibit 
spontaneous symmetry breaking [jTSj], the perturbative expansion breaks down and thus 
the necessity to implement resummation techniques. 

In this contribution we summarize the effects that the introduction of a finite pion 
chemical potential, associated to the finite pion density, has on the dispersion curve of 
pions at finite temperature. Starting from the linear sigma model, we use an effective 
Lagrangian [ [T4l ] obtained by integrating out the heavy sigma modes and compute, in 
a non perturbative fashion, the pion self-energy. We find that the dispersion curve is 
modified with respect to the vacuum in a way described by the introduction of an index 
of refraction larger than one, in addition to the thermal and density increase of the pion 
mass. Further details can be found in Ref. [Ol 



EFFECTIVE LAGRANGIAN 



The Lagrangian for the linear sigma model, including only meson degrees of freedom 
and with an explicit chiral symmetry breaking term, can be written as [[T^ 



where K and o are the pion and sigma fields, respectively, and the coupling is 



ml -ml 



When interested in a given approximation scheme to build the Green's functions of 
the theory at a given perturbative order, it is possible to exploit the relations that chiral 
symmetry imposes among them. These relations, better known as chiral Ward identities 
(ChWI), are a direct consequence of the fact that the divergence of the axial current 
may be used as an interpolating field for the pion [ITO]. In order to make sure that the 



approximation respects chiral symmetry, one needs to check that the modification of 
other Green's functions respect the corresponding ChWI. For example, two of the ChWI 
satisfied -order by order in perturbation theory- by the functions A^ (P) , (2) , T^j^ and 

n-^f' are 

04 

f^r^^{Q;0,P) = [A^'{Q)-A^\P)]5'J, (3) 

where momentum conservation at the vertices is implied, that is P^+ P2+ P^ = and 
Q + P = 0. The notation for the functional dependence of the vertices in Eqs. @ is such 
that the variables before and after the semicolon refer to the four-momenta of the sigma 
and pion fields, respectively [p^]. 

In Refs. [ [141 ] it has been shown that in the kinematical regime where the pion mo- 
mentum, the pion mass and the temperature are small compared to the sigma mass, the 
effective one-loop sigma propagator and one-sigma two-pion and four-pion vertices are 
given by 

= Q^-ml + 6X^fi^'{Q)- 
iFl^jiQ-P^^P^) = -2a^f^5'' [\-3X^y'{Q)] , (5) 

/r*;>«(;Pi,/'2,/'3,/'4) = 2a2{ 

X [l-3X^y\P^+P2)]5'j5^' 
+ [l-3X^y'{P^+P^)]5'''5j' 

+ [\-3X\y'iP,+p^)]d''dj'} , (6) 



where in the imaginary-time formalism of thermal field theory (TFT), the function J^'^ 
is obtained as the time-ordered analytical continuation to real energies of the function 
defined by 



^, ^ [ d^k 1 1 



(7) 



Here Q = ((U,q), K = ((OM,k) are Eucledian space four-vectors, namely = co^ + q'^, 
= (0^ + with CO = 2mnT and (0„ = 2nnT (m, n integers) being discrete boson 
frequencies, T is the temperature and (7 = |q|, = |k|. 

It is easy to check that Eqs. (Q), (|^) and (^ satisfy the Ward identities in Eq. (Q), this 
ensures that the approximation scheme adopted respects chiral symmetry. 

By using the above effective vertices and propagator, it is possible to construct the 
two-loop modification to the pion self-energy in the same kinematical regime with the 
result [0] 

Equation (§) reproduces the leading order result obtained from ChPT [Q]. Furthermore, 
we observe that Eq. can be formally obtained by means of the effective Lagrangian 

^='-{d,^)'-'-mlf-^_{f)\ (9) 

where a = 6{m\/2f]j.) and the factor 6 comes from considering the interaction of like- 
isospin pions in the vertex 

iYfi = -2i (^^^ (^5'^' 5^' + d'''5J' + . (10) 

Eqs. (|8|) and (^) mean that in the kinematical regime where the sigma mass is large 
compared to the pion mass, the momentum and the temperature, the linear sigma model 
Lagrangian reduces to a 0^ Lagrangian for effective like-isospin pions with an effective 
coupling given by a. In essence, the theory thus constructed and summarized by the 
effective Lagrangian in Eq. (|^) can be thought of a theory for the effective coupling a. 
By restricting ourselves to the above kinematical regime, we will proceed on working 
with the Lagrangian given by Eq. @. 



NON-PERTURBATIVE PION SELF-ENERGY 

It is well know that quantum field theories at finite temperature present certain subtleties 
such as the breakdown of the perturbative expansion [|17]]. This breakdown becomes 
manifest in two important cases: the appearance of infrared divergences in theories 
with massless degrees of freedom, and the compensation of powers of the coupling 



constant with powers of T for large temperatures. In both situations, the resummation 
of certain classes of diagrams represents an important improvement for the study of the 
physical properties of such theories. Even for cases where neither there is a massless 
degree of freedom, nor the temperature is extremely large, it is important to consider a 
resummation scheme, particularly for the case of theories with spontaneous symmetry 
breaking near critical behavior. 

In order to consider the above mentioned general situation and with the purpose of 
studying the pion dispersion curve at finite density and temperature, let us formally 
consider the Schwinger-Dyson equation for the pion self-energy 

^/^N f d^k 1 



2 {iTtfK^^ml+n 



'l'"2 



1 1 

K^ + ml + n{K^+K^-P)^ + ml + n' ^ 

where for internal lines, we make the substitution ico ico + /i. Notice that since the 
interaction Lagrangian contains only a quartic term, there is no need to dress the vertices 
in the above equation. 

Equation dTT] ) represents an integral equation for the function Tl(P), which, needless 
to say, is very difficult to be solved exactly. In order to find an approximate solution let 
us write 

n(/')~no+n(/') (i2) 

and consider tl{P) <^ TTq. As we will see, such assumption is justified given that in our 
approximation, fl ~ ^(a^). Keeping only the lowest order contribution in fl, we find 

^ _oc^ f d^k 1 



CX 2 {' d k Y d ^2 1 
1 1 

(14) 

Equation (|13]) represents a self consistent equation for the (momentum independent) 
constant TTq. This is the well known resummation for the superdaisy diagrams which 
constitute the dominant contribution in the large-A'^ expansion [|T3Jl of the Lagrangian in 
Eq. (|^). On the other hand, Eq. (|14|) represents a first approximation to the momentum 
dependent piece of n(/') . 



n(p) = 

X 




T(MeV) [I (MeV) 



FIGURE 1. ^mi + n^ as a function of (a) the temperature T for different values of the chemical 
potential ranging from /i = to /i = 130 MeV, from bottom to top, and as a function of {b) the chemical 
potential pL for different values of the temperature ranging from T = 50 to T = 150 MeV, from bottom to 
top. 



Figure 1 shows the behavior of the quantity y m| + Oq as a function of (a) the 
temperature T for different values of the chemical potential /i and {b) as a function 
of jU for different values of T. In both cases, y^'m^ + ITq grows monotonically with both 
T and ^. 

The pion dispersion relation is given as the solution to 

pI-[p^ + ml + IlQ + KelV{pQ,p)\ =0 (15) 

for positive Pq. In Eq.([T5|), fl'' denotes the retarded version of ll upon analytical contin- 
uation to Minkowski space. The left panel of Fig. 2 shows plots of p^ as a function of 
p for different values of /i and a temperature T = 120 MeV. As can be seen from this 
figure, YI{pq,p = 0) contributes to the increase of the pion mass. Also, for large p, the 
dispersion curves approach the light cone, always from within the causal region Pq> p^. 

We also notice that the dispersion curves can be parametrized by a function of the 
form 

PQ = ^n-^{T,^)p^+M^{T,^). (16) 

This is shown in the right panel of Fig. 2 where the dots represent the computed values 
of the dispersion relation and the continuous curves the fits. We thus observe that the 
presence of a finite chemical potential has two dramatic effects on the behavior of the 
pion dispersion curve compared to the case where only the temperature is considered: 
first, there is a significant increase in the pion mass and second, the index of refraction 
parameter n becomes larger than unity. 



SUMMARY AND CONCLUSIONS 



In this paper we have considered the effects of a finite pion density on the pion dis- 
persion curve at finite temperature. The finite density is described in terms of a finite 
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FIGURE 2. Left panel: as a function of p for ^ = 0, 120, 130 MeV, from bottom to top, and T = 120 
MeV. Right panel: Fit of the dispersion curve to the functional form = y^n^^{T,ii)p^ +M^{T,IJ.) for 
jJ. —0, 120, 130 MeV, from bottom to top, and T = 120 MeV. The dots are the computed values of the 
dispersion relation. 



pion chemical potential. We have argued that such description is important during the 
hadronic phase of a collision of heavy nuclei at high energies between chemical and 
thermal freeze-out when the strong pion-number changing processes have driven the 
pion number to a fixed value. 

In order to consider a general scenario that takes into account resummation effects, 
we have presented an approximate solution to the Schwinger-Dyson equation for the 
momentum-dependent pion self-energy, writing this as Tl{P) — IlQ + tl{P) and con- 
sidering tl{P) <^ Dq, which is justified given that in our approximation, tl ~ ff{a^), 
whereas the perturbative expansion of Hq starts at ff{a). 

The pion dispersion relation thus obtained at finite density and temperature deviates 
from the vacuum dispersion relation and can be described in terms of a density and 
temperature dependent mass and an index of refraction, also temperature, density as 
well as momentum dependent. This index of refraction is larger than unity for all values 
of the momentum. Similar results have been obtained in Ref. [ jT9[ ] using also a linear 
sigma model without the introduction of a finite chemical potential but rather as a 
consequence of the loss of Lorentz invariance when a particle travels in a medium at 
finite temperature. We note however that our result is more general than that of the 
former reference where the computation of the pion dispersion relation was carried from 
the onset in the weak coupling regime, that is to say, ^ 1, whereas our approach is 
valid for arbitrary values of i^M- 

Recall that one of the most salient features resulting from the analyses of collisions 
of heavy nuclei at SPS energies is the low mass dilepton spectra which shows an in- 
crease around the free p peak. Such behavior has been explained in terms of possible 
in-medium modifications of the p with either a shift toward lower values of its mass 
or an increase in its width []20|]. Given that the p meson couples strongly to the two 
pion channel, any possible in-medium modifications of the pion propagation properties 
should translate into induced modifications of the p meson properties and should there- 
fore be properly accounted for in a self-consistent manner in these kind of analyses. In 
particular, the increase of the pion mass with density modifies the phase space for the 



production and decay of p's in equilibrium and consequently it affects the production of 
e+ pairs. This matters will be the subject of a future work pT|]. 
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